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INTRODUCTION 

The condensation of water vapor in an air stream has the following 
consequences: 

1. Acquisition of heat (liberated heat of vaporization) 

2. Loss of mass on the part of the flowing gas (water vapor is 

converted to liquid) 

3. Change in the specific gas constants and of the ratio k of 

the specific heats (caused by change of gas composition) 

A discontinuous change of state is therefore connected with the 
condensation; schlieren photographs of supersonic flows in two- 
dimensional Laval nozzles show two intersecting oblique shock fronts 
that in the case of high humidities may merge near the point of inter- 
section into one normal shock front. The following discussion will 
deal with normal shock fronts only; it will be assumed that the 
velocity vector may be considered as being at right angles to the 
shock front (one-dimensional theory).- All flow properties directly 
ahead of the shock will be designated by the subscript 1, all those 
directly behind the shock by the subscript 2, all those referring 
to a stagnation condition by the subscript o, and all those applying 
in the narrowest cross section of the nozzle F* by an asterisk. 

[NACA comment: Particular attention is called to the fact that 

capital T denotes stagnation temperature only when used with the 
subscript o. Otherwise it denotes static temperature.] The sonic 
velocity will be denoted by a, and the flow velocity by w. 

The following equation exists for the sum of kinetic and caloric 
energy (enthalpy i = CpT) per unit mass of the flowing gas: 


*Analytische Behandlung des geraden Kondensationsstosses . 
Heeres-Versuchsstelle, Peenemunde, Archiv Nr. 66/72, Marz 30, 1942. 
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( 1 ) 


Directly "behind the shock the constant total energy i D has "been 
supplemented "by a quantity of heat energy per unit mass q as follows: 



i 0 + <1 


- I k + 1 a *i 2 
" 2 k - 1 a 


d 1 ) 


Because i c ' > i D , for the de s cription of the 'flow behind the shock 
a higher stagnation temp e rature" must be specified than that which 
characte fifed the flow before the shock . 

It will become apparent later that the values of p c , p D , and 
F* must also be changed. From equations (l) and (l 1 ) for a*, is 
obtained 

s Q > 1 . (2a) 

a* 

a * * 

as Q was thus defined as an abbreviation for another equa- 

tion may be written 

^1, Q2 (>1) (2b) 

x o J-o 

From equation (2b) it follows that with i 0 ' = i Q + q 


Q 2 = 1 + f- (3) 

1 o 

The quantity Q (inherently positive) is clearly determined by the 
quantity of heat of vaporization liberated; Q increases with that 
quantity. 

The loss in mass of flowing gas due to condensation, as well as 
the alteration of the gas constants and of k will be assumed to be 
so small that they may be disregarded. Other conditions being the 
same (for example, relative humidity = 100 percent), this assumption 
becomes more exact in proportion as T 0 is lower. If the a_r flows 
into the Laval nozzle from the free atmosphere with T 0 ar.nroxi- 
mately 273° K, the gravimetric proportion of saturated water vapor 3s 
2 to 3 percent (this proportion remains the same through all changes 
of state so long as these changes occur without gain or loss of water 
vapor ) . 
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Accordingly , the following equations apply: 

Pl v 1 = pg w 2 (4) 

Pq + Pi w x 2 = p 2 + p 2 w 2 2 (5) 

Pl Tp ' P2 T 2 W 

The equations (4), (5), and (6) express the constancy of the 
mass flow, the momentum, and the specific gas constants R. Together 
with equations (l) and ( 1 * ) , these equations constitute the basic 
equations of a simplified theory of the normal condensation shock. 

If the condensation occurs in a subsonic flow, the associated pressure 
change, because it is propagated at sonic speed, quickly spreads 
throughout the whole flew. The same pressure spread applies to all 
aspects of the change of state. No shock phenomenon can develop. 
Therefore, in the following discussion a fundamental assumption is 
that w^ is a supersonic velocity. 


CHANGE OF STATE AT SITE OF SHOCK 


A suitable combination of the five basic equations leads through 
a rather long mathematical process, which will be omitted here, to 
the following relation: 


p 

Wo W]_ w 2 — a*' 

W 1 W 1 w 2 “ a* 2 


(7) 


For a*’ 


a* (Q = 1) 


this becomes the Prandtl formula 


wq w 2 


= a* £ 


By introduction of the abbreviations 


W]_ Wp 

— rr = x n and ~~ 
a.* 1 a*' 


= x.- 


equation (8) is obtained 


w- 


Wi 


- Q 




( 8 ) 


Fx-om equation (7) is obtained 
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Q = - 


Xl x 2 - Q 


X 1 X 2 “ Q 


or 


2 1 + *1 
x 2 - x 2 ■ - Qx^" - + 1 “ 0 


(7a) 


From these equations, an equation to he used subsequently is obtained 

Qxi (1 + x 2 2 ) = x 2 (1 + xq 2 ) (7b) 


and by further conversion 


v 2 w 2 2 + a*' 2 
W 1 " wq 2 + a* 2 


(7c) 


which is another form of equation (7) 


The solution of equation (7a) defines x 2 as a function of Xq 
and Q. With the aid of this equation, simple expressions for p 2 /pq, 
T 2 /Tq, Ma 2 /Maq, ar.d p 2 /pq may be derived, [NACA comment: Ma, 

Mach number.] 

From equations (4) and (7) it follows that 

(9) 

Pi Q x 2 

If the well-known relation between T/T 0 and w/a* is applied 
to Tq/T 0 and T 2 /T 0 ' , by consideration of equation (2b) also, the 

following equation is obtained: 


t 2 

T, 


= 


, k-1 2 

1 ~ k?r x 2 


k-l 

k+2 


X 1 


( 10 ) 


Mag _ 'f? _ Y. 2 / T ±. 

Miq ~ w{ a 2 “ Wq '\j T 2 


Because furthermore, 
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it follows that 


Mag x 2 


Man 



k-1 2 

k+1 X 1 

k-1 2 

iSI x 2 


( 11 ) 


From equation (6) and with the aid of equations (9) and (10) 

2 


P2 

Pi 


Q 


x-, 1 


k-1 

k+1 


c 2 q _ 


2 2k 

x 2 “ 1 + x i - i^r ^i x 2 

___ _ 2 = , kjgL 2 

k+1 X 1 1 k+1 X 1 


k-1 


( 12 ) 


The expression on the right is derived from the middle expression 
by the application of equation (7b). 


The solution of equation (7a) is 

1 + Xq 2 ± / {l + Xq2)2 _ 4Q2 Xl 2 


Xg =, -- 


2QXq 


Because Xg, a ratio of velocities, must be a real number, from the 
previous equation an upper limit is obtained for Q; the lower limit 
is by definition Q = 1. It is evidently true chat 

isesL sr- <13) 

n 

Because in the supersonic region —■ > 1, accordingly the absolute 
value of the root is found to be 

(xq 2 - l) s? v /"(l + xq 2 ) 2 - 4Q 2 xq 2 ^ 0 

In order to decide the sign of the root, Xq being fixed and Q, 
being variable, the values Xg may assume for each of the two signs 
must first be determined. For the variation of xg with Q' it 
appears that 

&Xg _ ± Xg (1 + Xq 2 ) 

^ ‘ q7(i + 

where the absolute value of the root is in the denominator. With 
the minus sign, xg is a continually decreasing function of 3„, 
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whereby it follows from equation (13) that 

Xp = x 2 = 1 (14a) 

Correspondingly, with the plus sign (x 2 increasing with Q) it is 
found that 

< x 7 < 1 (14b) 

x i = * * 

In the first case, the air flow directly behind the shock is always 
moving at supersonic velocity; in the second case, always at subsonic 

1 + Xp 2 

velocity. Only for Q = ! ^ m _ x = — — is exactly sonic velocity 
found in both cases. 


The first possibility is supported by experimental findings. In 
the schlieren photographs, Mach waves or shock fronts both of which 
can arise only at supersonic velocities may be observed behind the 
shock or at least at a certain distance downstream. If subsonic 
velocity existed directly behind the shock, a subsequent constriction 
of the cross section of the flow filaments would have to occur (in 
order to produce a minimum cross section). This course of the flow is 
improbable because the shock always occurs at a point where the walls 
diverge downstream. Therefore the positive sign must be chosen for 
the root 


1 + Xp 2 + 7(1 + x-l 2 ) 2 - 4Q 2 Xp 2 

X 2 = 2Qx7 


(7d) 


and equation (14a) is obtained. The normal conden sation sh ock effects 
a change from one supers onic ve locity to anothe r s u per so n 1c velocity . 
The nomograph (fig. 1) shows trie relation between Q, Xp, and x 2 
as determined by equation (7d). In regard to the end. point of the Q 
scale, see the remarks at the end of this section. 

After the sign of the root has been determined, x 2 = v 2 /a*' is 
clearly a function of Xp and Q. The same relation is therefore true 
of the quotients in equations (8) to (12). No change is effected in 
this respect if the quotients are multiplied by x-j , Pi/0 o > Tp/T^, 
Map, or pp/p 0 because these quantities are also clearly functions 
of xp = wp/a*. Because the quantity xp is in turn clearly a func- 
tion of the Mach number Map immediately ahead of the shook and Q 
is clearly a function of q in accordance with equation (3), the 
following principle is obtained: 
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Pressure, density, temperature, Mach number, and velocity imme- 
diately, behind a normal condensation, shock are clearly determined by 
the Mach number at which the shock occurs and by the quantity of heat 
liberated . 

Thus both these quantities must be known if determinations of 
the conditions directly behind the shock are desired (assuming that 
P 0 , Po> To; and a* are known). Instead of the Mach number, the 
given quantity may, of course, be w^/a*, Pj_/p 0 ; Pi/p 0 ; T l/ T o; or 

F*/F because in the supersonic region all these quantities may be 
expressed as functions of each other. 

At what Mach number (or temperature) the supercooled vapor 
condenses, in what quantity the water is precipitated, and how much 
heat is accordingly supplied to the air stream, are questions that 
the resources of gas dynamics are Insufficient to answer; hence, 

Ma^ and q must be regarded here as variables that are independent 
of one another. However, in accordance with equation (13), there is 
for every value of x^ (or Ma-j_) a maximum for Q, and hence also 
for q. Therefore, If the condensation occurs at a certain Mach 
number Map, the liberation of more heat than a certain maximum 
d etermined by Map is impossible. From equations (13) and (3) 
is obtained 


, 2 . .2 

- (xq - 1) 

-0 ^ q i i 0 L 

4x-j_ 2 


(13a) 


For xq = 1, q = 0. Hence the normal condensation shock neve r 
occurs in the narrowest cross section of the nozzle ; supersonic 
velocity must first have been attained. 


With consideration of equations (13) and (14a), the limits may be 
determined between which, with a given value of x^ (that is, with a 
given site of the shock), the quotients in equations (8) to (12) must 
lie. From what has been said in connection with equation (14a), it 
is evident that X£ decreases as Q increases. In accord with equa- 
tion (8), the same relation is also true for wg/wq, whereby 



1 



(8a) 


In accordance with equation (4), the quotient p /p 
contrary manner 1 


acts in the 
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Pp 2x- 

1 £ — ^ 2 

1 + 


(9a) 


It is evident from equation (10) that Tg/Ti increases simultaneously 
with Q; that is 


To 2 

1 — < — £ — 
“ T x - k + 1 


(1 + x-, 2 ) 2 


4X! 




k-1 2 

k+r x i , 


(10a) 


The ratio 


Ma 2 v 2 


Mai wi 1/ T 2 

Wg /W]_ and T^/Tg decrease). 


decreases with increasing Q (because 


1 ^ 


Ma c 


- Ma 


1 


2= A 


/k+1 k -1 2 

W X") 


x l c 


1 

Ma-]_ 


(11a) 


Because p 2 /p^ an ^ Tg/T^ increase, pg/Py a T so must increase 
in accordance with equation (6). If Q = Q^ y , the root in equa- 
tion (7d) is eliminated and the equation becomes: 

2c Wx x 1 x 2 = 1 + x l 2 

From equation (12) is derived 


1 < — < 


p 2 L 1 + X 1 


Pl - k + 1 k-1 2 

1 “ k+1 X 1 


(12a) 


The limits in the left side of equations (8a) to (12a) are valid 
for Q = 1, that is, for the case in which no condensation shock is 
present at all. It is evident that: 

The normal condensation shock effects an increase in pressure, 
density, and temperature, whereas velocity and Mach number diminish . 

The ordinary normal compression shock, which occurs when a flow 
meets superior pressure, has these same characteristics. Nevertheless 
the compression shock must be regarded as something different in 
principle from the condensation shock, which is initiated by the intro 
duction of a quantity of heat. The difference is shown by the fact, 
among others, that the compression shock always produces subsonic 
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n 

velocity. Also in such a case the ratio wg/w-j = 1 /x-]_ decreases 

with .increasing xp, whereas in the case of the condensation shock 
it increases, as may be seen from the derivative that may be obtained 
from equation (7d): 

/ w 2 ^ / x 2 \ 

*\*lj „ H*l/ 2 Q(Qx 1 -x 2 ) 

“sr~ = Q 


~5xi 


A 


> 0 


. 2 x 2 ..2 2 

1 + ' X]_ ) — 4Q x 1 


because 


Q > 1 and 


xi > x 2 


The same statements apply for the absolute value of the velocity 
wg . In regard to the compression shock, this follows from the Prandtl 
equation wg = a*/x-^; in regard to the condensation shock, from the 
fact that if *2/ x l increases with xp, Xg and therefore Wg = a*'Xg 

acts in the opposite manner. 


must also increase . 


In the case in which 


The quotient pg/pq 

IYTT 


a* 


and 


Pi = Pi = T 1 


= 0, 


the 


k - 1 ” r '-L ~ “ > 

formula presented in equation (7) is inapplicable. For as soon as 
this extreme condition is approached, the concepts and therefore the 
equations of continuity that were used in the derivation of equa- 
tion (7) no longer apply. 

If this inapplication is disregarded for the moment, equation (7) 
is valid for all shocks occurring at finite nozzle cross sections 


( X 1 < A / 


k + 1' 


;x i < 


The condition expressed in equation (13) in conjunc- 


tion with equation (7) gives 

Q-max 


- + X-i 

X 1 J 


a function that increases with x-, if Xq > 1. If the quantity xq 
is allowed to approach its maximum value r— qy-p an (unattainable) 
absolute upper limit for Q, is obtained a3 follows 


^ma 


x max 


/k~ 


1.424 


1 
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To this value there corresponds an absolute maximum value for the 
quantity of heat liberated; from equation (13a) is obtained 


^max 


max 


k 2 - 


~ 1.027 i. 


Thus no more heat than this can be added by a normal condensation 
shock under any circumstances. If for example T 0 = 273° C, with 

c-n = 0.238 — the following is obtained: 

* gram degree 0 


Imax max 53 S6 . 7 


calor ie 

gram 


Of Xj_ . 


In figure 2, Qmax and in figure 3, (q/i D ) are shown as functions 


RATIOS p 0 '/p 0 , Pq'/Po, AHD F*’/F* 

Increase of Entropy 

It was pointed out in the Introduction that for the description 
of the flow behind the normal condensation shock a stagnation temper- 
ature T 0 ' > T 0 must be specified. In the following discussion, it 
will be shown, first by means of a special case, that F*, p Q , and 
p 0 also must be altered. 

It follows from equation (lla) that for Q = Q max the following 
equation applies 


= ( Ma 2)min = 1 

In this special case, exactly sonic velocity occurs directly 
behind the shock. Because this occurrence, however, is always present 
in the narrowest cross section of a flow, the shock cross section 
Fg > F* is therefore to be regarded as the minimum cross section with 
respect to the flow behind the shock. 

On the same basis, the relations must be written in this case 

P 2 = p*' and p 2 = p* 1 

From equations (9a) and (12a) it also follows, that if P]_/P 0 and 
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P]_/p 0 are expressed in termB of x-^ according to known procedure 
and if the expressions for P*/P Q and P*/P D are introduced 


p*' 





1 _ 

k-1 


1 + x n 


>* 4 


j . 
\k~l 


rrr x i) (1 + *i 2 ) v* 4 & 


Because a normal isentropic flow exists behind the shock (if a second 
condensation is excluded). 


1 



therefore 



(15a) 



(15b) 


The relations that have thus been found for Q = Q max , namely, 

Po' 4 P 0 > Po' ^ P 0 » and F *' > F *> wil1 be P roved for Q ! Qmax 
also and at the same time will be more precisely expressed. 

If Pq'/Po equals K and the calculations are started from the 
identity 
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by use of known gas -dynamic "formulas and equation (12) the following 
expression is obtained: 



In view of equation (7d), it is evident that the ratio p 0 ! /p 0 is a 
function of Q and Xq, that is, the ratio depends upon the quantity 
of heat q liberated by the condensation and upon the Mach number 
directly ahead of the shock . Thus the ratio can have the value of 1 
only for certain special values of these variables. In order to obtain 
some basis for the determination of the possible range of values of K, 
the behavior of K with the variation of Q and Xq will next be 
investigated . 

With regard to Q, (xq being held constant), it is sufficient to 
confine the investigation to the expression 


K = 1 - 


k - 


k 

k-1 


k + 


T x 2 2 ) ( x + x i 2 “ FTT ^i x 2 


Then 


3K _ 2k / _ k_q 
5Q k + 1 \ k + 


2k-l 

- 1 2\ k-1 j^ X 2 

t x 2 ) isr 


C 2 (1 + Xq 2 ) - XqQ (Xg 2 + l) 


(^ k - 1 v 2 

1 ~ nr x 2 


x l x 2 ^ ~ k + 1 A 2 


5 } 


By application of equation (7b), the quantity in the brackets becomes 
0 and therefore 


Bk 


< 0 


The ratio p 0 '/p Q decreases (at a given value of Wq/a*) with increasing 
Qi The maximum value of E is thus reached at the minimum value of 
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I\ 


Q = 1. In this case, a*' = a* and x-^xg = 1- Thus from equa- 
tion (16) is obtained 


*Wx ~ 


l -^ x 2 

2k 1 k+1 X 1 


X 1 ~ 


(17a) 


This formula is naturally the same as that applying to the ordinary 
normal compression shock for p Q '/p 0 . The maximum value that K ^ y 
can assume through the variation of x-^ is accordingly reached when 
x-j_ = 1 with K^y mn y = 1. This may also he mathematically shown 


dK max _ _ 2k 2-k 

dx, k + 1 • K max 


(x x 2 - 1) 
2 

.. 

k+1/ 


/k + 1 

Because l|x^S J p the derivative is negative (except at the 

limits of the range of values of x^ where it disappears). For 
x-] = 1 there is thus a maximum value of (namely, K^y = l) 

V k + 1 

a minimum value (namely, ^ = 0), which at 

the same time is the smallest conceivable value for K (because 
p 0 '/p 0 cannot become negative). Hence 


The ratio of stagnation pressure Po'/To lies between 0 and 1 and 


assumes the second value only if the condensation shock is nonexistent. 


B ecause in the deri vation of this theorem the sign of 
J(l + x-j^) 2 — 4QX]_ 2 plays no part, it is valid for both signs.) 
Accordingly it is permissible to describe K as a throttling factor 
as is customary in compression shock theory. 

When Q = Q^y, a minimum value of K is reached for the given 
value of x^ 

1 

? r ~ik-i 


^ain 


1 + X-|_ 
2 


k + 1 
2 


k - 1 
k + 1 


(17b) 
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This relation is the same as equation (15b). In this case, the 
velocity directly behind the shock wg is equal to the sonic veloc- 

1 + x-i 2 

ity a*' . The inequality 1 1 Qg — 2x~ corres P onds to tlle 


inequality > K ^ K^. 


Figure 4 shows the quantities K 1nay and K m j_ n as functions of 

x-^. When x-^ = quantities become equal to 0. However, 

in accordance with the remarks at the end of the last section, this 
result is only of theoretical interest. 

Now, conversely, Q rather than x-^ will be considered as given 
and fixed. It is found that 


(\ - l£=k v 2 

5K _ 2k \ k+1 1 . 

dx-^ k+1 


2-k 

k-l r 


\ k+1 


2k-l 

k-l 


SI i 1 - Irr *i 2 ) to (1+ ^i 2 )- (1+ x 2 2)] 


(1 - F7T x 2 2 ) ( x i 2 - 1 )(Q x 2 - x i j 


By application of equation (7b), the quantity within the brackets 
becomes 0. Furthermore, because from equations (8) and (8a) 


( Vr 


^2 ~ X 1 = X 1 W[~ 1 J <0 


for all values of x^_ > 1 


5k 

dxi 


< 0 


The ratio p 1 /p Q (at a given value of Q) decreases with increasing 
Mach number. 


Thus according to this reasoning also, the throttling factor 
assumes the theoretical minimum value K m in^ = 0 corresponding to 
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/k + 1 


( x iU 

with, the given value of 
tion (13) 


At the smallest value of x^ that is in agreement 
Q, K "becomes largest. According to equa- 


\ (k * "i) £ 


Q 


Because for x^ > 1 the expression on the left side decreases with 
decreasing x;j_, the smallest possible value of Q corresponds to the 
smallest possible value of x-^ denoted by W. Accordingly, W is 
determined by the equation 


1 + V £ 
2W 


= Q 


whereby 

W = Q ± Jq£ - 1 

1 + W^ 

If — gy — is substituted for Q, it becomes evident that (because 
W > l) in order to create the identity 

W = Q + /q 2 - 1 


if the root is understood to signify its absolute value. When this 
value is inserted in equation (16) (x 2 becoming equal to 1) the 
following is obtained: 


K mai Q) 


i + (q +7q 2 - i) 2 
2 


+ 1 


1 - 


k - 

k + 


Y (q + i ) 5 


1 

k-1 


(19) 


The right side of equation (19) is formally the sa me as that of equa- 
tion (17b) if x-^ is replaced by W = Q + Jo? - 1. The mathematical 
basis for this is that with a fixed value of x-j_ the quantities %a,x 
and x-j_ are so related by the same equation as the quantities Q 
and (x]_) mirL = ¥ with a fixed value of Q, that in both cases the 
expression for K in equation (16) undergoes the same simplification. 
Directly behind the shock, sonic velocity exists in both cases; this 
velocity thus occurs in every case wherein either the maximum quantity 
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of heat is liberated at a given Mach number, or the quantity of heat 
being given, the condensation takes place at the lowest possible Mach 

number. Figure 5 shows as a function of Q. 


Because it Was assumed (see Introduction) that the specific gas 
constant E is not altered by the shock, an equation may be written 


Po 


Po 


1 


T ' 

o 


From this equation and equation (2b) it follows that 

Po' _ _1_ V < 1 
P 0 Q 2 Po 


( 20 ) 


because it is true that Q > 1 and K < 1. The ratio of the stagna- 
tion densities cannot be greater than 1. 


Furthermore, because the loss of mass in the flow due to the 
precipitation of the condensed water is to be considered insignificant 
and disregarded, the following statement of continuity is valid: 

F *.p*i a *« - F *p*a* 

By consideration of equations (15a) and (2a), the following equation 
is obtained: 


F*' 


= F* 


Pq 1 
Po' * 


and because of equation (20) 


F*' 


F* Q 


Po_ 

Po' 


> F* 


( 21 ) 


The second equality sign is valid only when Q = 1 and p 0 '/P 0 = 
that is, when there is no shock at all. For the description of the 
flow behind the normal condensation shock, a larger minimum cross 
section must be specified . 

Being a discontinuous, that is an irreversible, process, the 
condensation shock involves an increase in entropy. If the entropy 
per unit of mass of the flowing gas ahead of the shock is written as 


NACA TM No . 1174 


17 


a = c p In T q - H In p D 

and that behind the shock is written as 

s' = Cp In T 0 1 — B In p D ' 

with the aid of equation (2b) the following equation is obtained: 

As = s' — s = Cp In — B In K (22) 

In accordance with equation (16) the following statement applies: 

The increase of entropy due to the normal condensation shock is a 
function of the Mach number at which the shock occurs and of the 
quantity of heat liberated . In the case of an ordinary normal com- 
pression shock, the first member on the right side of equation (22) 
is lacking; consequently, with a given throttling factor, the increase 
of entropy due to a compression shock is smaller than that due to a 
condensation shock. 


FLOW PROPERTIES AT EXIT CEOSS SECTION OF NOZZLE 

The flow properties in the exit cross section Fg of a wind- 
tunnel nozzle are the same as those in the whole test section of the 
tunnel and are therefore of especial interest. All flow properties 
referring to Fg will be denoted by the subscript E; all flow 
properties that have undergone alteration due to a condensation shock 
will be denoted by a prime (except that those directly behind the 
shock will have the subscript 2 as previously). Behind the shock a 
normal isentropic flow exists to which the known gas-dynamic laws can 
be applied. These laws connect the conditions at any selected cross 
section F downstream of the shock with the conditions directly 
behind the shock; these conditions are therefore functions of 
Xg = Wg/a* and Q. Hence Fg, Xg, and Q are determinative of 
the conditions at the exit cross section Fg . If two of these 
quantities are held constant, the conditions at the exit cross section 
depend only upon the third quantity. 

As previously shown, the quantity P 0 /p 0 ' increases both with 
fixed Q, and increasing xg (a transfer of the shock downstream) 
and with fixed Xg (fixed shock site) and increasing Q. Therefore, 
if Fg and Q are held constant (or Fg and Xg), in accordance 
with equation (21), F*' and therefore F*'yFg also increases with 

increasing Xg (or with increasing Q) . If, however, F*'/Fg 
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increases, then in the region of supersonic flows (which is what is 
actually found behind the shock) pg ' /p Q ' , pg'/p 0 , and T e '/T 0 * 

increase, while w-g'/a*’ and the Mach number Mag ’ decrease. 

With a fixed exit cross section Fj; and a given value of F*, 
the flow properties in the nozzle are determined. 

If Q also is regarded as fixed in accordance with equations (2a) 
and (2b) neither a*' nor T 0 ' are altered. Thus Wg' decreases 
and Tg 1 increases with increasing Xg. From the continuity rule, 
Pg'w-g' = constant, it follows furthermore that pg' increases and 
from the equation of state ITi'/Pe'^e' = con8 ‘ tan ' fc that pg 1 also 
increases. Because Mag increases with Xg and q remains consts.nt 
with Q (compare equation (3)), the conclusion is made that: 

If in a nozzle the Mach number at which a certain specified 
quan tity of heat o f condensation is supplied to the flow is increas ed 
(that is, if the shock occurs f urther downstream), at the exit cross 
section the Mach number and the absolute value of the velocity are 
less but the pressure, density, and temperature are greater . 

If in addition to Fg, the Mach number Mag is considered (that 
is, xg also) as fixed, then according to equtation (2b) T 0 ' increases 

with increasing Q and Tg' thus increases all the more. However, 
a*' also increases while Wg'/a*’ decreases. Hence no general state- 
ment concerning wg' is possible. This impossibility is also evident 
from another line of reasoning. Namely, if the exit cross section and 
the location of the shock coincide Wg' = w 2 ; according to equa- 
tion (7d), w 2 (which equals x 2 Qa*) is a decreasing function of Q. 

If on the contrary the exit cross section is assumed as infinitely large 

V = a *Q VlHHr 


where wg ! - is an increasing function of Q. In accordance with the 
continuity rule, the behavior of the density with increasing Q thus 
also depends upon the size of Fg, that is (at a given F*), upon 
the Mach number for which the nozzle is constructed. In order to 
investigate the variation of pg* as a function of Q, the calcula- 
tions may again be started from the equation of state by writing 


Pe’V - P*V 


Pe* V 

p* v 
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From Figp^'v-g 1 = F*p*a* is derived 


Re ' F* a*_ 

P* _ F E w e ' 

and with T 0 ' = T 0 (compare equation (2b)) is obtained 

- , T . _ 0 *t 11 q 2 a^_ V 
?E 1 E - P i 0 p ^ w ' T ' 

E E o 

Because according to equation (2a), Qa* = a*', if T e '/T 0 ' is also 
expressed in the known manner in terms of w E , / a *' 


Pe'T e ' - p*T 0 - 


Q 1-^(2. 

F* ^ k+1 \a* 


'w™ ' \ 2 


The equation of state thus assumes the form 


= constant 


f ^E * ’ 

\a*' < 


Differentiation with respect to Q yields 


a,* 

= constant — 


W E k-1 / W E 

1 k+1 Va*' 




_ 0 1 + 3 Eri 

L_ k+1 \a*V; 


( v jlL) 

\a*V 


Because wp'/a*' is a decreasing function of Q, 5 1 ^ — J ■ 

Spp* \PQ/ 

-g— > 0 where pg.’ increases with Q. 


V\ 


< 0 and 


The over-all result is as follows: 


If in a nozzle, the site of the shock being fixed, the quantity 
of heat of condensation liberated increases, the Mach number at the 
exit cross section decreases while the pressure and the temperature 
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Increase. The behavior of the densit y and of the a bsol ute value o f 
the velocity depends upon the Mach number for -w hi ch the' wind tunnel 
is built . ~~ " ~~~~ “ ' 

If a normal condensation shock occurs In a nozzl e , the same 
statements apply . [KAOA comment: It is believed that this sentence, 

completely redundant in the present form, was intended to read: "If 

an oblique condensation shock occurs...”.] As a consequence of the 
shock, q increases from q = 0 to a value q > 0. 

In the absence of a shock, a Mach number equal to the reduced 
Mach number that would exist at the exit cross section if a shock had 
occurred will occur at a point in the nozzle upstream of the exit 
cross section. The quotients p/p D , p/p 0 , T/T Q , an(i w/a*, which 
are clearly functions of the Mach number, must at that point have the 
same value as the corresponding exit cross section quotients pg'/p 0 ’, 


Pe7 p o'> 

that 

VAo'. 

and 

w E 7 a *' 

if a shock occurs. Hence 

it follows 

V _ V 

p = Po 

.< 1 , 51 

’ p 

_ P °' 
= Po 

T-o 1 

■<1, -f 

T n ' 0 w-p 1 

• - fy - 9 > 1 ’ and v- ■ 

Rif* 

5*-= 9>1 


That is, 

If in a nozzle, as a consequence of a normal condensation shock, 
the Mach number at the exit cross section is reduced, at the exi t 
cross section the pressure and the density are lower and the tempera- 
ture and the absolute velocity are higher than w oul d cor re spond to 
the same Mach number in the absence of the shoo k (this Mach number 
would then occur at a cross section upstream of the exit cross section) . 


If the last two of the three quantitites Fg, x^_, and Q are 
held constant, that is, if the site of the shock and the quantity of 
heat of condensation liberated are regarded as fixed given values, 
then with increasing Fg pressure, density, and temperature naturally 
decrease, while the Mach number and the absolute velocity increase, for 
downstream from the shock a normal isentropic flow exists. Answering 
the question of how the quotients pg'/pg, pg'/pg, t e 7 T E; Ma E ! /Mag, 
and wg ' /wg behave as functions of Fg is not equally easy. 


By consideration of 
ance with equation (2a), 


the last of these quotients first, in accord - 
the following equation is applicable: 


V 

V E 


= Q 


V 

a* 1 


wg 

a* 
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from which hy formal development 



Introduction of the known general relation between F*/F and w/a*, 
a process of computation, which will here he omitted, leads to 



According to equation (21), F*'/F e > F*/F e and therefore 
w E '/ a * 1 < w E /a* and the bracketed expression is thereby greater 
than 0. The quotient w E ' /w E thus increases with F E ; it assumes 
its least value if Fj; coincides with the shock front cross section, 

and its greatest value if F T — »«. In the first case, , 

w E ‘ w E /vTTT We W1 

and m the second case ^7- = ~ = /w . By consideration of 

equation (8), the following equation is obtained: 


Q 



g Q 


(23) 


The gas velocity at the exit cross section can increase up to Q 
times the original velocity.. 

The expression on the left side of equation (23) is equivalent 
to Wg/w-j^ < 1 (compare equation (8 a)); that on the right side is 
Q > 1* Equation (23) explains the fact already mentioned that, 
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assuming F* to be given, upon the occurrence of a normal condensa- 
tion shock the behavior of the velocity at the exit cross section 
depends upon the magnitude of this shock. The velocity decreases if 
the exit cross section decreases and increases if the exit cross 
section increases. This relation is evident. The shock effects so 
great an increase of pressure, density, and temperature that the 
velocity is first decreased. Because the flow contains a greater 
total amount of energy after the shock, with a sufficiently great 
expansion, that is, with a sufficiently marked decrease in the pro- 
portion of caloric energy, the velocity must finally become greater 
than it could have been without the shock. 


In order to determine the exit cross section Fg at which the 
velocity would remain precisely the same, the calculation may be 
started from the expression 



According to equation (21), F*'/F* =• Q/K. By inserting the relation 

Wg 1 = w-g = Wg, the following equation is obtained: 

/*E\ 2 _ k + 1 Q 2k - Q 2 !^" 1 
Va */ " k - 1 Q 2k _ K k-1 


The quantity K is a function only of the fixed selected values 

and Q. Corresponding to Wg/a* in accordance with the known relation 

a certain value of F*/Fg exists from which Fg may be calculated . 

/wg\ 2 k + 1 

Because Q > 1 and K < 1, < k — p that is, F E has a finite 

value . 


From the continuity expression pg'Wg' 
is obtained 


PgWg and equation (23) 


1 *1 Re' 1 

Q x 2 - Pe ' Q 


(24) 


In the cross section Fg, jog' = pgj upstream of that point, 
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Cg' > pg, and downstream pg' < pg. In spite of the compressing 
effect of the condensation shocks the density at the exit cross section 
can he less than it would have been in the absence of the ' shock. 


As a consequence of the energy equation 

V * 1 + S 

the following equation may he written 

Wip * W-jji 

~ + V " ~ 2 ~ + % + q 

where with i = CpT it follows that 


w- 2 /’w® 1 2 


/ w r 

W 


- ! I = Cp Tg 


1 - 


L J£_ 

% 


+ <1 


Because Wg’/vg continuously increases with increasing Fg, v/% 
must continuously decrease. The least value of this quotient is 
reached when Wg’/wg = Q. For this case, from the preceding equation, 
If q is expressed in terms of Q in accordance with equation (3), 
the following equation is obtained: 

\ (a 2 - i> + »p % (55- - i) - io - D 


Because 

Wg 2 

2 + °p % = ^o 

it follows from the previous equation that 



The maximum value for Tg'/Tg is T 2 /T^ (condensation shock occurring 
at the exit cross section). Therefore the following equation is 
obtained (compare equation (10)): 


1 - 


1 - 


k-1 2 
k+T *2 
k-1 2 
k+1 X 1 


= T- 


V 


E 


> Q 


(25) 
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As may "be seen, the statement is confirmed that, as the result 
of a normal condensation shock (Q 4 l)-, the temperature at the exit 
cross section is raised. 

For the ratio of the Mach numbers, the following relation is 
obtained : 


Ma E* _ fE V E ' pE~ 

Ma E Wg a E ' w E y T E ’ 

Thus the ratio increases with F E and as F E — » co the ratio approaches 
the value 1, in accordance with equations (23) and (25). The least 
value is that of Ma 2 /Ma]_. Accordingly, by use of equation (ll) the 
following equation is obtained: 


X 1 


h ~ 2 

,:L k+1 X 1 


Ma E ' 


k-1 2 a* Maw - 

1 - hi *2 * 


< 1 


(26) 


This equation corroborates the fact that the normal condensation shock 
reduces the Mach number at the exit cross section. 


The equation of state yields 

PE* _ PE^_ . 

PE Pe t e 

The pressure ratio is, like pg'/pg and V /T e , a decreasing func- 
tion of F e . By use of equations (12), (24), and (25) the following 
equation is obtained: 


y 2k 

1 + *1 - gx 0*i x 2 Pe' 


k~l 


> 




(27) 


It is evident that the normal condensation shock increases the pres- 
sure at the exit cross section. 


In equations (23) to (27), x 2 may in all cases be clearly 
expressed in terns of x^ and Q, in accordance with equation (7d). 
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The results previously derived are valid without change in 
respect to any nozzle cross section downstream of the shock F that 
may he chosen to substitute for Fg. Therefore the results may be 
summarized as follows: With increasing distance downstream of the 

site of a condensation shock, the ratios p'/p, p'/p/ and T ! /T 

decrease while Ma’/Ma and w'/w increase. The numerical values 
of the flow properties constituting the numerators and the values 
of those constituting the denominators vary in the same manner. 
Therefore it may be concluded that: 

Downstream of a normal condensation shock, the pressure, the 
density, and the temperature decrease more sharply with increasing 
nozzle cross section and the Mach number and the velocity increase 
more sharply than in the absence of the shock. 


Translation by Edward S. Shafer, 
National Advisory Committee 
for Aeronautics. 
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